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Abstract 

We study the 3d Ising universality class using the functional renormalisation group. 
With the help of background fields and a derivative expansion up to fourth order we 
compute the leading index, the subleading symmetric and anti-symmetric corrections 
to scaling, the anomalous dimension, the scaling solution, and the eigenperturbations 
at criticality. We also study the cross-correlations of scaling exponents, and their 
dependence on dimensionality. We find a very good numerical convergence of the 
derivative expansion, also in comparison with earlier findings. Evaluating the data 
from all functional renormalisation group studies to date, we estimate the systematic 
error which is found to be small and in good agreement with findings from Monte 
Carlo simulations, e-expansion techniques, and resummed perturbation theory. 

1. Introduction 

Continuous phase transitions of numerous systems in statistical and particle physics 
belong to the Ising universality class, characterised by the short range nature of the 
interaction, a scalar order parameter, and the dimension. The absence of a physical 
length scale at the phase transition implies scale invariance. Many fluids, magnets, or 
particle physics models thus share the same behaviour at criticality described by universal 
numbers such as the scaling exponent for the correlation length u, its subleading correction 
u, and the anomalous dimension of the order parameter at criticality t]. Further critical 
exponents for eg. the specific heat a, the spontaneous magnetisation 0, the magnetic 
susceptibility 7, and the magnetization at criticality as a function of the external field 5, 
are linked to v and r\ by scaling relations 

The computation of universal indices - in a quantum field theoretical or statistical 
physics setting - has become a benchmark test for perturbative and non-perturbative 
methods in field theory. Accurate predictions for exponents, scaling functions or ampli- 
tude ratios are available based the renormalisation group, resummations of perturbation 
theory, and numerical simulations on the lattice (see {2J for an overview). An important 
continuum method in the above is the functional renormalisation group [3|, |4[ , based on the 
infinitesimal integrating-out of momentum modes from a path integral representation of 
the theory with the help of a Wilsonian momentum cutoff [5] . By construction, functional 
flow equations continuously interpolate between the microscopic theory at short distances 



and the full quantum effective theory at large distances. Powerful optimisation techniques 
are available to maximize the phy sics content in given approximation, and to minimize cut- 



off artifacts along the flow [6H10] . A particular strength of the functional renormalisation 



group is its flexibility, allowing for perturbative as well as non-perturbative approxima- 
tions even in the presence of strong correlations or couplings 1CH1 
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Fixed point studies for Ising-like theories have been performed within various reali 



sations of the functional renormalisation group including Polchinski's formulation 
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imation 
expansion 



Wetterich's equation [12J, exact background field flows [16Hl9|, the proper-time approx 
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, and discretized (hierarchical) transformations 25 



26 . The derivative 



28[ and variations thereof |l2|, |29| are the expansion schemes of choice 



based on a small anomalous dimension. Scalin g be haviour in more complicated theo 



ries, eg. thermal field theory |30[, gauge theories [31n33| and gravity [34J, can equally be 
accessed using thermal or gauge- covariant derivative and vertex expansions. 

In this paper, we study the Wilson-Fisher fixed point in three dimensions within a 
background field formulation. In the past, background field methods have mostly been 
employed for gauge theories and gravity, where they allow for a gauge invariant implemen- 



tation of the cutoff , [35I, [ifl], also offering new expansions schemes The motivation 
for using this technique for non-gauge systems is two-fold. First of all, the presence of 
a background field allows for a re-organisation of the flow equation. While this is of no 
relevance for the full flow, it does make a difference once approximations are invoked. 
In particular, derivative expansions of standard and background field flows are different. 
This allows for complementary measurements of universal scaling exponents. Secondly, 
background field flows have provided very good numerical results to lower orders in the 
derivative expansion. Therefore it is important to understand whether this pattern carries 
over to higher order. 

In addition, we discuss the convergence of the derivative expansion and provide an 
estimate for systematic uncertainties. Error estimates are obtained by probing the de- 



pendence on the shape of the Wilsonian cutoff function 38|, |39] - which vanishes for the 
physical theory and hence should become small with increasing order in the expansion 
- and by checking the numerical convergence of successive orders, which we extend up 
to fourth order. We estimate the systematic error by comparing different projections of 
the Wilson-Fisher fixed point onto the flow equation, using a weighted average over the 
available data. We find a coherent picture and very good agreement with the mean values 
and error estimates from Monte Carlo and perturbative studies. We also evaluate the 
cross-dependences of scaling exponents, and find an interesting link between the expected 
error in an observable and its sensitivity to tiny variations of the dimensionality. 

The format of the paper is as follows. We recall the basic set-up (Sec. 2) and our 
approximations (Sec. 3), followed by a discussion of results (Sec. 4) and their optimisation 
(Sec. 5). Two sections deal with the dependence on dimensionality (Sec. 6) and the 
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cross-correlation of exponents (Sec. 7). We evaluate the convergence of the derivative 
expansion (Sec. 8) as well as systematic uncertainties (Sec. 9), and close with a brief 
discussion (Sec. 10). 



2. Renormalisation group 

Wilson's renormalisation group is based on the integrating-out of momentum de- 
grees of freedom from a path integral representation of the theory. Modern, func- 
tional, implementations of this idea employ an infrared momentum cutoff term ASk = 
\ J ip(q) Rk{q) <f(—q) for the propagating modes <p(q), added to the Schwinger functional 
with classical action 5* and external current J, 

In Z k [J] = In j [dip] exp f-S -AS k + J if ■ J^j . (1) 

The cutoff function Rk{q) can be viewed as a momentum- dependent mass term with k 
denoting the RG momentum scale. It obeys Rk(q 2 ) — > for k 2 /q 2 — > to ensure that the 
large momentum modes q 2 ^k 2 can propagate freely, and Rk(q 2 ) > for q 2 /k 2 — > which 
ensures that the low momentum modes q 2 ^ k 2 are suppressed in the functional integral. 
This makes Rk an infrared cutoff. The change of (OQ) with the RG scale k (t = In k) reads 
d t Z k = —(d t ASk)j- In terms of the effective action r fc [0] = supj(— In Z k [J] + • J) + ASk 
it is given by Wetterich's flow equation ji| 

d t r k [cf>] = ^Tr-^ d t R k , (2) 

an exact, functional differential equation which links the scale-dependence of Tk[4>] with its 
second functional derivative 1^ [0] = 5 f^f^ and (the scale-dependence of) the regulator 
function Rk- The trace denotes a momentum integration, and <fi = (<p)j denotes the 
expectation value of the field (p at fixed external current J. By construction, the flow (T2]) 
interpolates between an initial microscopic action I\ w S at k = A and the full quantum 
effective action at k — 0. 



Next, we discuss background field flows following |16Hl9l| where a non-propagating 
background field is introduced into the effective action Tk[4>] — > ^k[4>, 4>] by coupling the 
fluctuation field (0 — 0) to the regulator and the external current. For the derivation of 
the flow, the background field acts as a spectator, and we obtain ()2]) with the replacement 
ri 2 ' > [0] — > 6 • The background field dependence of r fe [0, 0] is governed by 

-=F k 0, = -Tr— T — (3) 

W 2 rf[0,0] + R k 50 

and vanishes in the infrared limit k — > where R — > 0. Subsequently, the background 
field will be identified with the physical mean = 0, leading to a background field flow 
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for an effective action T k [(p} = T k [(f>, <$>]• This technique is standard practice in the study 
of gauge theories and gravity leading to gauge-invariant flows within the background field 
method. In gauge and non-gauge theories, this procedure can simplify the evaluation of 
the operator trace in (]2]), which makes it attractive for our purposes 18|, |24j. The key 
difference between standard and background fields stems from the fact that the presence 
of the background field, at an intermediate stage of the computation, corresponds to a 
re-organisation the flow. This aspect is exploited below. 



To be specific, we introduce the background field by substituting q 2 — > T 



(2), 



in the regulator function R k (q 2 ] 
For some class of ^-functions 



(other choices such as R k — > R k {4>] can be used as well). 
18], the flow (P2]) takes a very convenient form, 



d t T k = Tr 



T^/rn + k 2 



o{d t v 



(2h 

k ) 



(4) 



Here, m G [1, oo] parametrises a remaining freedom in the choice for the cutoff function 
which we fix later. The term oc d t T k originates from the implicit ^-dependence introduced 
in Rk via T k , and reflects the re-organisation of the flow through background fields. The 
term d t T k ' on the r.h.s. of (jl]) can be replaced through a series - starting off with the 
leading term in fll]) and functional derivatives thereof - by making repeated use of (jl]). 
Closed forms for <9tIV or the flow for a few relevant couplings are available under certain 
approximations 0, Eii 32]. Below, we need the flow for several field-dependent functions 



and therefore limit ourselves to the leading term [121, [18 . 

The first term on the r.h.s. of (HI), or linear combinations for various m, is equivalent 
to Liao's proper time flow equation [13, Esi 
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d t T, = - 



iTr / ™(d t f k )exp 
Jo s 



, 6 2 r k \ 



(5) 



originally derived from a proper time regularization of the one-loop effective action. 
Eqs. (|4]) and ([5]) are linked via f k = / PT (sA 2 ) - fp<r(sk 2 ) with /pt(z) = T(m,x)/T(m). 
The flow (|1]) can equally be obtained from generalized Callan-Symanzik flows [l7| without 
the necessity for background fields. The flow equation (HI) in the approximation flS]) has 
oreviously been used for studies of phase transitio ns Il6l . l2lH24|, tunneling phenomena 



40| | . spontaneous and chiral symmetry breaking 4ll443|. gravity 4J], and a general proof 



of convexity 19]. Here, we will use it to analyse the infra-red scaling at the Wilson-Fisher 
fixed point to fourth order in the derivative expansion. 



3. Approximations 

In this section, we detail our ansatz for the effective action for a real scalar field based 
on the derivative expansion and the relevant renormalisation group equations. With up 
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to fourth order derivative operators, the effective action reads 



d D x 



V k {4>) + \Z k (<f>) d^d^ + W k (<j>) (<9 2 0) 2 . (6) 



The ansatz ([6]) should capture the relevant infrared physics provided the anomalous dimen- 
sion of the fields stay small. Note that the derivative expansion has no small parameter 
directly associated with it, because the integrand of ([2j) receives dominant contributions 



for q 2 /k 2 ^ 1 28] . Hence the numerical convergence has to be checked a posteriori. Good 
numerical convergence is known for appropriate momentum cutoffs 

The three functions V, Z and W in our ansatz dSJ) are symmetric under reflection in 
field space <p <H- —<p. In principle, there are three independent tensor structures available 
to fourth order in the derivative expansion, 

WM (<9 2 0) 2 , H k {</>) d,<pd^(d 2 <p) , J k (<j>) {d^d^f (7) 

with J (H) symmetric (anti-symmetric) under reflection in field space. In the free theory 
limit, the operators ([7]) scale identically. At an interacting fixed point this degeneracy 
is lifted, and the higher derivative operators contribute with different strengths to the 
flow (j3J). We expect that the term ~ W is the most relevant one, for reasons detailed 
in Sect. 4. Therefore we neglect H and J. Then the initial conditions for the flow at 
momentum scale k = A are 

y A (0) = im 2 2 + iA A 4 , Za(0) = 1, W A (<f ) ) = 0. (8) 

For k < A, higher-order couplings are switched on due to the renormalisation group 
running (j4j), and the functions V, Z and W develop a non-trivial field dependence. The 
Wilson-Fisher scaling solution for k — > corresponds to critical initial conditions m\ c and 
Aa, c - The renormalisation group equations for the functions V, Z and W are obtained 
by inserting into (j3j) and expanding the exponential by making use of the Baker- 
Campbell-Hausdorff formula. The partial differential equations for V, Z and W are of 
the form 

a ' x r d i i (9M e " A ' Kx <9) 

where X — V, Z or W, and Aq = V" + Z p 2 + 2 W p A , with primes on functions denoting 
derivatives w.r.t. the fields. The equations (Q encode the central physics of our setup. 
The kernels Kx encode the interactions amongst the operators in the ansatz ([6]) under 
the renormalisation group. We have Ky = 1. The kernels Kz {Kw) are polynomials in 
the loop momentum variable p up to order p u (p 20 ) with coefficient functions depending 
polynomially on V, Z, W and their derivatives, and the proper-time integration parameter 
s. The expressions are very long and not given explicitly. 

For a fixed point study, it is convenient to introduce dimensionless variables and a spe- 
cific cutoff. We mainly use the parameter 1/m = which is equivalent to the step function 
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approximation 


x 


v"(0) 


v"(x Q ) 


z(xq) 


w(0) 


LPA 


1.899 


-0.297 


0.672 


1 





o(d 2 ) 


1.889 


-0.266 


0.601 


1.047 





o{d A ) 


1.888 


-0.267 


0.609 


1.050 


-1.26 10~ 4 



Table 1: Potential minimum xq, curvature, and other reference values of the scaling solution. 

fpT(y) = OJl -y) with y = sZk 2 , achieved as /pt(J/) = lim m ^oo Y(m, my)/Y{m) 
(see also [45]). For large rj^ 3> k 2 , the flow then becomes exponentially suppressed 
oc exp(— Y^ /k 2 ), rather than algebraically. In consequence, amplitude expansions con- 



vergence more rapidly 28|. The remaining s-integration in is performed analytically. 
We introduce 

* = In(Jfe/A), x = k 1 - 012 - 7112 , p = p/k, a = A k- 2+ ^ 
v(x) = k~ D V(<P), z{x) = k r ' Z{<f>) , w(x) = k 2+r > W{<t>) , (10) 

where we have rescaled dimensionful variables in units of k; rj = —dt In Zk denotes the 
anomalous dimension of the field. It is understood that v, z and w are functions of t and 
x. Below, we denote derivatives w.r.t. x as eg. d x v = v' . In the parametrisation (fTUj) . the 
explicit /c-dependence of the differential equations (Q is factored into the variables. We 
finally obtain 

d t Y + D Y Y-D x xY'= f A-$ e -(f*l*) k y with Y = {v',z,w} (11) 

J (2vr)^ 

The terms on the left-hand sides display the canonical and anomalous scaling of the fields 
D x = [4>] and the variables Dy, with 

£>* = §(£ -2 + 77), D« = \{D + 2-rj), D g = - V , D w = -(2 + V ) . (12) 

We note that the scaling dimensions D x and D v > (D z and D w ) are positive (negative) 
for 77 > and D > 2. The terms on the right-hand sides parametrise the non-trivial 
interactions induced by ([6]) under the renormalisation group. The integral kernels K Y are 
related to the kernels Kx\ s =i in © via the relations K v i = —d x (a /z), K z = Kzk D+v , 
and K w = K w k D+2+ri . 

4. Results 

In this section, we analyse the physics of ( ITTjl at the Wilson-Fisher fixed point, which, 
in D = 3 dimensions, corresponds to the unique non-trivial solution Y*(x) 7^ of 

d t Y* = 0. (13) 
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approximation 


n 


V 


LO 


L0 5 


A 


A 5 


LPA 





0.6260 


0.762 


2.163 


0.477 


1.354 


o(d 2 ) 


0.0330 


0.6244 


0.852 


2.459 


0.532 


1.535 


o{d A ) 


0.0313 


0.6247 


0.865 


2.563 


0.540 


1.601 



Table 2: Anomalous dimension, leading and sub-leading scaling exponents, and the Wegner 
corrections A = vlo and A5 = voj§ for different orders in the derivative expansion (see text). 

In the limit of large x ^> 1, the r.h.s. of (fill is exponentially suppressed and the fixed 
point solution is dominated by the scaling of the fields and variables, 

Y*(x) oc x Dy/d * . (14) 

Consequently, the solutions z*(x) and w*(x) vanish asymptotically because D z and D w 
are < 0. The algebraic suppression is the more pronounced the larger —Dy. For v', we 
find a rising behaviour for large x because D v > > 0. For small x ^ 1, the interaction terms 
become relevant, and the complexity of the equations makes it necessary to use numerical 
methods. Here, we solve ( ITTjl with ( 1T31) for v[(x), z*(x) and w*(x) without making any 
further expansions such as eg. polynomial expansions. 

The results are displayed in Fig. [IJ Including the wave function renormalisation z*(x), 
the first derivative of the potential v^(x) changes only mildly from the local potential 
approximation result. The further inclusion of w*(x) leaves v#(x) practically unchanged, 
whereas the wave function renormalisation z*(x) increases mildly, though only for larger 
x. We note that w*(x) is very small and negative for small x, enhancing its impact for 
smaller x. Some characteristic values of the scaling solution are given in Tab. [U Including 
second (fourth) order operators, the vacuum expectation value changes approximately by 
1% (0.1%), the curvature v" by 10% (1.5%), and the wave function renormalisation by 
5% (0.5%). 

We point out that the quantitative relevancy of operators in the effective action corre- 
lates with their scaling dimension. This is already visible from the results to lower orders 
in the derivative expansion. We have D u > > > D z > D w , which materializes at the 
fixed point ( fT3j) as variations in the scaling solutions u'^, z* and of order 1, 10 _1 , and 
10~ 3 , respectively (see Fig. CQ. Quantitatively, this can be understood as follows. The 
exponential suppression of terms on the right-hand side of ( fTTj) for large field variable 
x implies that the large-field behaviour of operators is solely determined by their mass 
dimension, see Q14p . The transition from small-field behaviour to large-field asymptotics 
is exponentially strong, thereby restricting the impact of higher derivative operators to 
the small field regime. The same observation applies for the variation of the fixed point 
solution and for the eigenperturbations at criticality under successive extensions from 
V(<f>) -> V(<j>), Z(<p) V(<P), Z{(f>), W{(j>). 
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Figure 1: Wilson-Fisher fixed point in D = 3 for (a) the first derivative of the potential v'^{x), 

(b) the deviation of the wave function renormalisation from its classical value 10 (z*(x) — 1), and 

(c) the four-derivative operator 10 3 u;*(a;). Coding: local potential approximation (black dotted 
curve), 2nd order derivative expansion (red dashed curves), 4th order derivative expansion (blue 
continuous curves). 



Next, we comment on the approximation ([6]). A full 0{d 4 ) order calculation in the 
derivative expansion of the effective action requires the inclusion of the terms H and J, 
see flTJ). Close to the free field theory limit, the three terms in Q scale identically, but 
this degeneracy is lifted at a non-trivial fixed point solution. The mass dimensions of 
W,H and J in (J7|) are different, and increasingly negative, eg. D w = [W] = —(2 + 77), 
Dh = [H] = — ^(5 + 377) and Dj = [J] = — (3 + 2rj) in D = 3 dimensions. Hence 
> D w > D h > Dj and, consequently, the scaling solutions h*(x) and j*(x) will be 
suppressed compared to w*(x), see (JHj). Therefore we expect that the impact of H and 
J on scaling exponents is subleading, analogous to the pattern observed to lower orders 
in the derivative expansion. We also note that the suppression, in general, will depend 
quantitatively on the regularisation. The suppression is exponential for the background 
field flow used here, and hence stronger than the power-law suppression observed for 



standard flows in a derivative expansion, see |46M48l| . 

Small deviations from the fixed point <&y,\(x) = Y(x) — Y*(x) are classified according 
to their universal scaling exponents A. In the vicinity of the fixed point the eigenpertur- 
bations obey the eigenvalue equation 

d t ® Y ,\ = A $y, A • (15) 

We solve (fl5j) using ffTTj) and the fixed point solution ([TBI to find the leading and sub- 
leading eigenvalues as well as the eigensolutions. The leading eigenperturbations $(x) are 




Figure 2: Eigenperturbations $>(x) with eigenvalue v at the Wilson-Fisher fixed point for (a) 
the potential v^(x), (b) z*(x) (rescaled by a factor 10), and for (c) w*(x) (rescaled by a factor 
10 3 ). Coding as in Fig. [TJ 
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Figure 3: Eigenperturbations <3?(x)with eigenvalue to at the Wilson-Fisher fixed point for (a) the 
potential v'*(x), (b) z*(x) (rescaled by a factor 10), and for (c) w*(x) (rescaled by a factor 10 3 ). 
Coding as in Fig. [TJ 
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info 


7 


IS 


V 


a 




5 


UJ 


a) world average 


1.2372(5) 


0.6301(4) 


0.0364(5) 


0.110(1) 


0.3265(3) 


4.789(2) 


0.84(4) 


b) this work 


1.2298 


0.6247 


0.0313 


0.1259 


0.3221 


4.818 


0.865 


c) impl. opt. 


1.233 


0.627 


0.034 


0.119 


0.324 


4.803 


0.839 



Table 3: Comparison of a) the world average of theoretical predictions [2J, with our results b) 
(|16p and c) f)19[) . also using scaling and hyper-scaling relations. 



symmetric under x -H- —x. The eigenvalues obey Ao < < Ai < A2 < • • • with X = —1/is 
and Ai = u in the statistical physics literature. The eigenperturbations $(x) which are 
antisymmetric under under x <H> — x have eigenvalues < Ai < A2 < ■ • • , and the smallest 



eigenvalue is denoted as Ai = co^ in the literature; see [49] for a determination of C05 in the 
local potential approximation. Our results for the eigenperturbations are given in Figs. [2] 
and|3j The inclusion of z*(x) and w*(x) changes the eigenperturbations with eigenvalue 
v (iS) only mildly from the local potential approximation. For the scaling exponents and 
the anomalous dimension, we find 

v= 0.6247, 77 = 0.0313, to = 0.865 , u 5 = 2.563 . (16) 

The numerical precision for u (ujk) is of the order 0.1% (1%). Comparing with lower 
orders in the derivative expansion, Tab. [2J we conclude that the derivative expansion of 
the background field flow displays a very good numerical convergence. 

In Tab. EK) and b), we compare the 'world average' of scaling exponents in three 
dimensions as compiled in [2J with our findings ( jl~6{) . Most exponents agree on the 
percent level and below, in particular the indices is, 7, 5 and /3 which are predominantly 
sensitive to the field-dependence of vertices at vanishing momentum. The anomalous 
dimension rj, and the exponents a and u are subleading and more sensitive to the 
momentum structure of propagators and vertices. Consequently, their precision is lower. 
Interestingly, the exponent u already agrees with the world average within 3%. The 
exponent a and the anomalous dimension rj only agree within 15% with the best values 
quoted in the literature. The same pattern persists the comparison with recent high 
accuracy results from Monte Carlo simulations is MC = 0.63002(10), u MC = 0.832(6), 
Vmc — 0.3627(10) (5o[. Here, the indices is, A = isoj, uj and r\ agree to within 0.8%, 3%, 4% 
and 13%, respectively. 

5. Optimisation 

In this section, we discuss the optimisation of results based on a background field flow. 
It is well-known that physical observables within an approximation of the functional flow 
can depend on the shape of the momentum cutoff function R and its parameters, eg. m. 
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The reason for this is that the cutoff function, a non-trivial function of momenta, couples 
to all operators in the effective action. Neglecting some operators means that some cutoff- 
dependent back-coupling in (|2j) is missing. Formally, within some systematic expansion 
of the flow equation, one obtains the exponents as a series 



phys = V(0)(R) + V(1){R) + V{2){R) + ■■■ , (17) 



where the contribution from every single order ur n -\ (R) may depend on the cutoff function 
R, and only the full physical result will be independent thereof. Optimisation is based 
on the observation that the convergence of ()17|) . and similarly for other observables, is 
improved through optimised choices for R, ie. the parameter m. 

To evaluate the -R-dependence of our results - and also as a consistency check for our 
numerical codes - we have re-calculated the exponents to leading and second order in the 



derivative expansion using (j4j) for other values of the cutoff parameter m |22H24j|. The 
output is displayed in Tab. H] and, as expected, the results fully agree with earlier findings. 
We add the following observations. To leading order in the derivative expansion, z/(m) is 
monotonous, covering the physical value v av « 0.63. This is different from the pattern 
observed using the standard flow in the same approximation (9)]. The index uj(m) is equally 
monotonous, approaching the world average value o; av ~ 0.84 from below, u(m) < w av . 
Together with the improved convergence of the amplitude expansion, this justifies the use 
of 1/m = 0, also favoured by a minimum sensitivity condition. We note that the physical 
value for v is matched at 1/m « 0.11. 

To second order in the derivative expansion, all three observables u, u and rj are 
monotonous functions of m. The ranges covered include the physical value in all three 
cases. A principle of minimum sensitivity is not applicable. For v (go), the relative change 
from leading to second order is minimal atl/m = 0(l/m~ 0.4). We can use our results 
to second order in the derivative expansion to identify the value of the cutoff parameter 
m = m av for which rj - the least well-determined index - matches best with the prediction 
from the 'world average' or Monte Carlo studies. We find l/m av ~ 0.08, which is very 
close to the best match 1/m « 0.11 found to leading order, and 

v = 0.626, 77 = 0.036, u = 0.823 . (18) 

The difference between (1181) and Monte-Carlo results for u, u and A reduces to 0.7%, 
1.6% and 2.3%, respectively. The improved agreement shows that the scaling exponents 
display the correct cross-dependences, also supporting small values for the inverse cut- 
off parameter 1/m. The fit l/m av ps 0.08 comes out slightly larger than 1/m = 0, a 
consequence of the anomalous dimension being underestimated in the latter case. The 
cross-correlation amongst u, u and rj as functions of the cutoff parameter m is similar to 



the strong cross-correlation observed earlier in the local potential approximation 26 



11 



m 


2 


2.5 


3 


4 


DC 


1 1 w J 


0.065 


0.056 


0.051 


0.045 


0.033 


z^Ilpa 

v \o(d 2 ) 


0.660 
0.632 


0.650 
0.631 


0.644 
0.630 


0.638 
0.629 


0.626 
0.624 


wlpa 

w \o(d 2 ) 


0.628 
0.677 


0.656 
0.702 


0.674 
0.725 


0.698 
0.756 


0.762 
0.852 



Table 4: Variation of scaling exponents with cutoff shape parameter m; background field flow 
to leading and second order in the derivative expansion. 



Using the above, we obtain an improved estimate to fourth order in the derivative 
expansion by combining (1161) with the m-dependence of indices for small 1/m from Tab. H] 
at 1/m ps 0.08. Note that at large m, the difference m^| (a 2 ) — m ^\o(d 4 ) m the m- 
dependence of exponents with X = u, u or rj between the 2nd order and 4th order results 
are small. With this approximation, we arrive at 

u = 0.627, 7] = 0.034, u = 0.839 . (19) 



The difference between ( fT6j) and ( 1T9|) serves as a measure for the variation of indices 
within the stable domain of RG flows to this order in the approximation. In Tab. Efc), 
we compile our findings ( TT9|) and compare with the world average. The agreement with 
Tab. [3|l) becomes significantly enhanced. For the indices 7, u, 77, a, (3, 5, A = u v 
and uj we find an accuracy of 3%, 0.5%, 7%, 8%, 0.3%, 0.2% and 0.6%, respectively. The 
same quality persists the comparison with recent high accuracy results from Monte Carlo 
simulations 50|. Here, the indices u, A = uu,u and rj agree to within 0.5%, 0.4%, 0.9% 



and 6%, respectively, which is a clear improvement over the results at 1/m = (see 
Sec. 4). 

In summary, optimisation of background field flows consistently favours small values 
for the parameter 1/m. The excellent agreement of indices around 1/m ps 0.08 indicates 
that the operators retained in our approximation display the physically expected 
cross-correlations. This non-trivial result lends additional support for the present set-up 
and the internal consistency of the underlying approximations. 



6. Variation with dimensionality 

In this section, we consider the Wilson-Fisher fixed point away from three dimensions. 
It is a useful consistency check to understand the global D- dependence of our findings, and 
their interpolation between the known results in two and four dimensions. Furthermore, 
probing the local D-dependence by perturbing the system with -j^|£) =3 provides insights 
into the structural stability of our set-up. We note that this information is also of interest 
for systems of finite size or finite volume, where variations of the (spatial) dimensions are 
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sensitive to the variations with linear system size. 

With increasing dimensionality D > 3, the scaling exponents approach mean field 
values at D = 4 with r\±& = 0, = \ and = 0. For decreasing D < 3, further higher- 
order critical points become accessible whenever n »s D /{D — 2) becomes an integer, with 
n = 3 corresponding to the Wilson-Fisher fixed point. In D — 2, scaling exponents and 
the anomalous dimension take the known values v 2 d = 1, ^2d — 1 and 772^ = \- 

We have computed u and rj in the vicinity of three dimensions, see Tab. EJ Below 
D ^,2.7, the identification of the scaling solution becomes numerically more demanding. 
This should be related to the appearance of a competing scaling solution, which becomes 
available when D »s 8/3. From the data, we find 

du do; dn , . 

_ = _0.18, ^ = -0.65, ^=-0.08. (20) 

for the first derivatives at D — 3. Note that the derivatives would read —0.25, —0.5 and 
—0.125 for a simple linear interpolation between the analytically known results at D = 2 
and 4. Fitting the data points for v and 77 (u) with a cubic (quadratic) polynomial in 
D, and extrapolating we find z^|d=4 ~ 0.49, uj\d=a ~ 0.02 and t)\d=a ~ —0.02. In the 
opposite limit, extrapolation leads to ^ |d=2 ~ 0.92, oj\d=2 ~ 1-3, and r]\ D=2 ~ 0.2. These 
estimates are fully consistent with the expected behaviour, and the slight deviations at 
the endpoints serve as (rough) indicator for the underlying error. The extrapolated result 
for r] is smaller than the exact one 77 = j, and the fourth order result is slightly smaller 
than the second order result. This suggests that our result slightly under-estimates the 
value for rj at D = 3, though a definite conclusion would require more data points for rj in 
2 < D < 3. For a study of 77 and v in 1 < D < 4 dimensions using an optimised standard 



flow to second order in the derivative expansion, we refer to [51 . 

Next, we estimate the relative variation of scaling exponents. Suppose we are interested 
in a physical observable X . The relative variation of X with D around the dimensionality 
of interest serves as an indicator for the stability in the observable X. A low stability 
indicates that an observable will depend more strongly on the approximation, and vice 
versa. With this in mind, we write dX/dD = —Ax X for X = u, u, rj. Using f )20|) . we find 

A v — 0.28(1) , A- = 0.85(1), A, = 2.5(1) (21) 

in three dimensions. The smallness of A v explains the high accuracy achieved for v 
already to low orders in a derivative expansion. In turn, the large value of A v explains the 
stronger sensitivity of 77 on the approximation. Furthermore, the pattern A v < A u < A v 
suggests that the expected accuracy in u should be better than the one in 77, and worse 
than the one in v. This is in accord with the pattern observed in our results, see Sect. 4, 
and with the earlier functional RG results discussed below (see Sect. 8 and 9). 
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0.762204 


0.0330 


0.0313 
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0.8899 


0.0519 
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2.7 


0.688 


0.949 


0.0637 
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Table 5: Variation of u, u and rj with dimensionality D to leading and second order in the 
derivative expansion (see text). 



7. Cross-correlations 

Cross-correlations amongst scaling exponents provide insights into the finer structure 
of the theory, and into the inner working of the approximation in place. Within the local 



potential approximation, cross-correlations are strong |26|, and only weakly dependent on 
the cutoff Rk, in particular for optimised flows {(J. A similar cross-correlation has been 
observed based on hierarchical RG transformations, thereby providing a link between the 
cutoff (Rk) dependence of the continuum RG and finite step size effects in discrete versions 
thereof 2a . 



Here, we are interested in the correlations to higher order in the derivative expansion. 
To set the stage, we perform a linear interpolation for the derivatives based on the known 
results at D = 2 and D = 4. We find 

£ = 2, £ = 0.5, £ = 0.25. (22) 
au au auj 

Within our functional RG set-up, we access the cross- correlation of exponents by keeping 
the regulator fixed, and by exploiting that ( 1201) represent full variations with D. Since 
rj(D) is monotonous in D, at least in the region of interest (see Tab. |5]), we invert rj(D) 
into D(rj) to obtain the functions 1/(77) = u {D(r})) ana w(rj) = oj(D(rj)) which encode the 
cross-correlation of scaling exponents. In three dimensions, their first derivatives read 

^ = 3.63, ? = 0.45, = 0.124. (23) 
au au au 

Note that ^f 1 ^ = 1 to within 0.03%, which is smaller than the error in (1231). 
Comparing ( 1231) with the linear approximation ( 1221) . we find that 4^ is roughly of the 
same size, while ^ (^) is roughly twice (half) as big as the linear approximation. Our 
result (123]) compares well with the estimate ^j| e - C x P . = 0.59 obtained from a modified 
epsilon expansion {52)]. We note that ^jIfrg < 7p\im. < j^le-exp.- The double- logarithmic 
derivatives follow from (121j) in an obvious manner, eg. = A^jA v ^ leading to the 
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estimates 44^ ~ 3 , ~ 9 and 4r^ ~ 3, consistent with the sensitivity observed in 

d m v ' din/ d In OJ ' J 

our results. 

8. Convergence 

Despite the small anomalous dimension, the Wilson-Fisher fixed point corresponds to 
a non-trivially interacting theory and is therefore intrinsically non-perturbative. While 
little is known about the absolute convergence of systematic approximations to fl2j) in the 
non-perturbative regime, the numerical convergence of expansions can be accessed order 



by order [28j. In this section, we discuss the convergence of the derivative expansion (see 
Tab. by comparing results for u, r\ and u from different realizations of the functional 
renormalisation group including the standard flow (J2J), background field flows (jlj), and the 
Wilson-Polchinski flow (see [28] for an earlier overview). We have omitted data points 
which are not based on an (at least partly) optimised choice for the momentum cutoff, 



eg. sharp cut-off results [531 ] . 

To leading order in the derivative expansion, the full cutoff dependence of v(R) is 
known within the standard flow 0, |26|, within the Wilson-Polchinski flow 15] where the 



result is unique, and, partly, within background field flows [2l|-|24|. For the standard flow, 



the best result is given in b) [9] , achieved for suitably optimised regulators. High-accuracy 



expressions for the exponents are given in 54J and are in full agreement with findings from 



the Wilson-Polchinski flow (l5j in c). The background field flow covers a larger range of 
values for z/(m), the smallest one given in a). Comparing a) with b), we note that the 
leading index v (subleading index oS) is slightly (significantly) closer to the physical result 
in the setup a). 

For the 0(d 2 ) approximation, we report the exp onents from the standard flow based 
on an optimised algebraic (power-law) cutoffs |271j in d), a standard exponential cutoff 



55l . |56| in e), an optimised exponential cutoff [47] in f), and a flat optimised cutoff [47 
in g). Note that algebraic (power-law) cutoffs of leads to slowly converging flows 
within the derivative expansion j^j, which is already visible within the local potential 
approximation j^]. The comparatively large estimate for 77 in e) is a consequence thereof. 
Below, we will retain e) for a conservative error estimate. Comparing d)-g) with h), we 
note that the indices v and u differ only slightly amongst the different implementations. In 
contrast, the anomalous dimension t] varies more strongly, about ±25%. In the standard 
flow, the anomalous dimension stays above 4%, whereas the background field flow leads 
to a result below 4%, closer to the physical value. 

Results to second order in the derivative expansion are also available within Polchinski's 
formulation of the renormalisation group 0, Et| . The Wilson-Polchinski flow is linked to 
(EI) by a Legendre transform, implying that derivative expansions are inequivalent beyond 



the trivial order. A significant cutoff dependence, in particular for 77, is observed 57- 
59], which calls for a stability-based optimisation of the cutoff js), [To|, 15]. A prediction 
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1) 


a 4 , bf 


0.0313 


0.6247 
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this work 
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<9 4 , bf, implicit 


0.034 
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0.839 
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0.039 
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Table 6: Comparison of results from the functional renormalisation group within various ap- 
proximations (see text). Local potential approximation (d°): a) background field flow (bf); b) 
standard flow (st); c) Wilson-Polchinski flow (WP). Derivative expansion to second order (<9 2 ): 
d) - g) standard flow (various cutoffs); h) background field flow. Derivative expansion to second 
order with matching of the anomalous dimension: i) Wilson-Polchinski flow; j) background field 
flow. Derivative expansion to fourth order (<9 4 ): k) standard flow; 1) background field flow; m) 
background field flow with implicit optimisation. Mixed approximation retaining momentum- 
and field dependences (mixed): n) standard flow. 



for Ising exponents is achieved at 0(d 2 ) by tuning the cutoff to the desired value for 77, 
and using a minimum sensitivity condition to identify the remaining 



say rj « 0.038 [58 
exponents 



see also 57, 59 



58 , summarized in 



This leads to v « 0.625 and u « 0.77 
Tab. ©). The predictions for v and u are in the expected range of values, showing that 
the Wilson-Polchinski flow displays the correct cross-correlation of scaling exponents. It 
will be interesting to see whether a fourth-order computation stabilises the result. For 
comparison, we have added in Tab. [6j) our result ( fl~8l) from the background field flow to 
second order in the derivative expansion, where 77 has been matched to the world average 
and Monte-Carlo result. The Wilson-Polchinski and background field estimates agree very 
well for is, and differ by less than 8% for the exponent 00. The background field value 
is much closer to the expected value. Note that this procedure is not applicable for the 
standard flow to second order, because the anomalous dimension stays above 4% for all 
cutoffs and cannot be matched to the physical value. 
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Beyond 0(d 2 ) in the derivative expansion, we cite the fourth order computation by 



48| in k), which is compared with our result (fT6l) in 1), the optimised background field 
result (Tl9|) in m) and the 'mixed' analysis of |60j in n), which is optimised using the 
principle of minimum sensitivity {si, |6l|, 62 1. The approach n) retains momentum- and 
field-dependences in the ansatz for the effective action, amended by approximations on 



the level of the flow; see |29[ for technical details. The results for v in all approaches 
are very close to the world average v &N = 0.6301(4). The value for u from background 
field flows 1) and m) are closest to the world average w av = 0.84(4). All values for r\ are 
now below 4%. Still, a slight variation of rj remains visible which makes the anomalous 
dimension the least well-determined obervable in Tab. [6j We note that the prediction for 
rj based on k) and n) are equally close to the world average rj av = 0.0362(4), approaching 
it from opposite sides. This is interesting because n) should have a better access to the 
momentum dependence of propagators. We suspect that the approximations on the level 



of the flow exercised in |29[ are responsible for this pattern. The rj- values from background 
field flows approach the physical value from below, with m) being closest to the expected 
value. 

The mean values based on all data points in Tab. [6] are (^)fRG = 0.630 and (w)fRG = 
0.790. For the anomalous dimension, we find (?7)fRG = 0.0312(0.0397), depending on 
whether we retain (suppress) the LPA data points ?7lpa = 0. (We come back to a detailed 
discussion of mean values and systematic errors in Sect. 9.) 

We use the numerical convergence of the derivative expansion for a crude error estimate. 
For the standard flow with order-by-order optimised exponential cutoff function Rk(q 2 ) oc 
a q 2 1 (exp q 2 /k 2 — l) we compare the LPA result r\ = and v = 0.6506 jfj with higher orders 
in the derivative expansion Tab. [6}:) and k). This leads to v « 0.637 ±2% and rj ~ 0.0387 
±15%. Retaining only the two best values for v improves the error estimate, v « 0.630 
±0.3%. The relative change Au/u reads 3.5 x 10 -2 (6.3 x 10~ 3 ) at second (fourth) order in 
the derivative expansion. For the background field flow with cutoff m — > oo, we compare 
Tab. Ek), h) and 1), leading to v « 0.625 ±0.4% and 77 = 0.0322 ± 2%. Hence, in the 
approximation (T4|), the numerical convergence of background field flows is slightly faster. 

We conclude that the derivative expansion of the functional renormalisation group, 
together with suitably optimised regulators, shows a very good numerical convergence 
up to high order for both standard and background field flows. 



9. Systematic errors 

Estimating systematic errors is common practice in eg. lattice approaches and re- 
summations of perturbation theory. Here we discuss how analogous estimates can be 
achieved for the functional renormalisation group, where physical observables are ob- 
tained by projecting the full flow in 'theory space' - the infinite dimensional space of 
operators parametrizing the effective action - onto a sub-set thereof. This step im- 
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plies an approximation and is a potential source for systematic errors. The flexibility 
of the formalism, however, allows for many different projections. Then the quantitative 
comparison of different projections gives access to the systematic uncertainty. 

We recall that, in general, approximations to the flow equation fl2]) enter via operators 
neglected in the effective action approximations on the level of the flow dtT^, and 
the choice for the momentum cutoff Rk- These aspects are partly intertwined, to the 
least because a momentum cutoff introduces a non-trivial momentum structure into the 
flow. In general, the operator content is central. A similar importance should be given 
to approximations on the level of the flow d t Tk, which feed back into the determination 
of scaling exponents. The regulator is crucial for the stability and convergence of the RG 
flow jsj. Within given approximations for Tk and dtT^, the regulator can be optimised to 
maximise the physics content in the flow, and to minimize cutoff artefacts. Uncertainties 
due to the boundary condition for the effective action are irrelevant for fixed point solu- 
tions. We conclude that systematic errors should only be derived from 'cutoff-optimised' 
results to eliminate cutoff artifacts 

Next we employ this reasoning to the data collected in Tab. EJ A first estimate for the 
systematic error is obtained by taking a weighted average over representative entries for 
each projection method (standard flow, Wilson-Polchinski flow, background field flow), 
disregarding further details of the approximations. Common to the data points is that 
the underlying regulators are, at least partially, optimised We first consider the data 
points Tab. |6k), b), f), h), k), 1), n) to obtain 

* = 0.631!*™ »? = 0.03618:881, ^ = 0.783^- (24) 

For rj, we only took the data with rj ^ into account. The mean values ()24|) change by 
less than 0.1%(1.5%) for u, r) (u) had we included the data points i), j) and m) based 
on some additional input. Hence ( T24"l) represents an average with equal weight for the 
different implementations of the functional flow. Note that the width of the error bars, 
roughly a standard deviation, are set by the least advanced approximations. 

An improved estimate is obtained by retaining only the most advanced results in 
Tab. O ie. k), 1) and n), all of which are based on a similar operator content, supported 
by a partial cutoff optimisation but differ in the projection technique. We recall that 
in k) a standard full fourth-order derivative expansion is used, together with a polynomial 
expansion in the fields (48| ; in n) a mixed approximation is employed retaining momentum- 
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and field-dependences, but neglecting loop momenta of certain vertex functions [60 
1) a background field flow is used within a fourth-order derivative expansion and without 
polynomial expansion in the fields, but neglecting higher order flow terms and subleading 
fourth-order derivative operators in the action (this work). The qualitative differences in 
the approximation make sure that the computations project in different manners onto the 
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resummed PT 

e-expansion 
world average 
Monte Carlo 



functional RGs 0.034(5) 0.630(5) 0.82(4) 



r\ v ijj ref. / year 

0.0335(25) 0.6304(13) 0.799(11) [52] ( i 9 98) 

0.0360(50) 0.6290(25) 0.814(18) [52] ( i998) 

0.0364(5) 0.6301(4) 0.84(4) [2] (2000) 

0.03627(10) 0.63002(10) 0.832(6) [5CT] (2010) 



this work 



Table 7: Comparison of results from the functional renormalisation group with resummed per- 
turbation theory, Monte-Carlo simulations, e-expansion, and a world average. 

Wilson-Fisher fixed point, thereby probing the systematic error. We find 

i/ = 0.630lffi, 17 = 0.03418:28, u = 0.823^. (25) 

Note that we omit the data set m) from this estimate to achieve a conservative error bar 
and an equal weight between projection methods. From ( 124]) to (|25|) the error bars are 
reduced by at least a factor of two. The mean values for z/, u and rj are shifted by 0.2%, 
5% and 6%, respectively. The shift in the mean values from (124"]) to ( 125]) is of a similar 
size as the estimated error in ([25]). 



In Tab. [7J the combined functional RG results (125]) are compared with the e-expansion, 
resummed perturbation theory, Monte Carlo simulations, and a world average of theo- 
retical predictions. It shows that the functional RG results agree very well with results 
from other methods within systematic errors and on the level of the mean values. The 
results are also compatible with recent experimental results, eg. 77 = 0.041 ± 0.005 and 
v = 0.632 ± 0.002 [63(, with experimental errors slightly larger than those from theory 
(see 64|, |65[ for overviews). Expected errors from the functional RG are presently about 
an order of magnitude larger than those from eg. numerical simulations, and more data 
and extended approximations are required to further reduce the systematic uncertainty. 
In particular, the value for u in ( 125]) is presently only based on two data points. Here, 
it would be useful to know the value from the standard flow at fourth's order in the 
derivative expansion to improve the error estimate in Tab. [7J Natural candidates for 
further data points are approximations of the Wilson- Polchinski equation beyond second 
order in the derivative expansion, or approximations with an improved access to the full 
momentum structure of propagator and vertices. 

10. Discussion 

The computation of universal scaling exponents is an important testing ground for 
methods in quantum field theory and statistical physics. We have obtained new results for 
the indices u, u, and rj of the 3d Ising universality class using functional renormalisation 
group methods within a background field formulation. Our analysis complements earlier 
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studies without background fields. Our findings to second (fT8j) and fourth order (fT6|) . (fT9l) 
in the derivative expansion agree very well with other theoretical studies. A fast numerical 
convergence of the derivative expansion is found, confirming similar observations to lower 
order in the expansion. The indices also display the correct cross-dependences. This 
non-trivial result lends further support to the underlying approximations. 

We have also studied the cross-correlations of exponents, and their sensitivity to tiny 
variations of the dimensionality. The latter correlates with the expected error of exponents 
within the derivative expansion. As a result ( l2Tj) . the index v shows a weak, the subleading 
index u a moderate, and the anomalous dimension a strong dependence on dimensional 
variations. We conclude that the achievable precision in these observables follows the 
same pattern, as confirmed by the data ( 1251) . 

The flexibility of the functional renormalisation group allows for different projections 
onto the Wilson-Fisher fixed point. We have exploited this freedom to estimate the 
systematic uncertainty of scaling exponents using all available data. The resulting mean 
values and error estimates (125)) agree very well with results from resummed perturbation 
theory and lattice simulations. More work and further data points are required to reduce 
the error bars, which are similar to those from experiment, but larger than those from 
recent numerical simulations. Natural candidates for further data points are eg. Wilson- 
Polchinski flows to fourth order in the derivative expansion, and approximation schemes 
with an improved access to the momentum structure of propagators and vertices. 

In addition, we have analysed the convergence of the derivative expansion, comparing 
data from standard flows, background field flows, and the Wilson-Polchinski flow. 
Background field flows lead systematically to smaller values for rj, and the derivative 
expansion converges very fast. Standard flows provide narrower bounds on exponents, 
while the derivative expansion shows a slightly slower rate of convergence. For the 
Wilson-Polchinski flow, structural arguments suggest that approximations beyond the 
leading order are more sensitive to the cutoff. Still, good results are available to second 
order, provided rj is matched. It will thus be interesting to extend these studies beyond 
the Ising universality class. 
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